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ABSTRACT 

We report the analysis of a highly magnetised neutron star in the Large Magellanic 
Cloud (LMC). The high mass X-ray binary pulsar Swift J045106. 8-694803 has been 
observed with Swift X-ray telescope (XRT) in 2008, The Rossi X-ray Timing Explorer 
(RXTE) in 2011 and the X-ray Multi-Mirror Mission - Newton (XMM-Newton) in 
2012. The change in spin period over these four years indicates a spin-up rate of 
—5.01 ± 0.06 s/yr, amongst the highest observed for an accreting pulsar. This spin-up 
rate can be accounted for using Ghosh and Lamb’s (1979) accretion theory assuming 
it has a magnetic field of (1.2±{j y)xl0 14 Gauss. This is over the quantum critical 
field value. There are very few accreting pulsars with such high surface magnetic fields 
and this is the first of which to be discovered in the LMC. The large spin- up rate is 
consistent with Swift Burst Alert Telescope (BAT) observations which show that Swift 
J045106. 8-694803 has had a consistently high X-ray luminosity for at least five years. 
Optical spectra have been used to classify the optical counterpart of Swift J045106.8- 
694803 as a B0-1 III-V star and a possible orbital period of 21.631±0.005 days has 
been found from MACHO optical photometry. 
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1 INTRODUCTION 

High mass X-ray binaries (HMXB) are binary systems com- 
posed of either a neutron star, white dwarf or black hole and 
an optical companion, either a supergiant star - in the case 
of supergiant X-ray binaries (SGXB) - or a dwarf, subgiant 
or giant OBe star - in the case of Be/X-ray binaries (BeXB). 
HMXB can be detected by the X-ray, optical and infrared 
emission they produce. 

X-rays are produced as matter is transferred from the 
optical star to its denser companion. It is accelerated as 
it moves into the gravitational potential well produced by 
the denser star and - if the denser star is a neutron star 
- it is then channelled by the star’s magnetic field lines to 
the magnetic poles. The accreted matter then rapidly de- 
celerates when it reaches the surface and potential energy 
is converted to heat which energises the plasma, producing 
X-ray hot spots. These can appear as ‘beams’ if they peri- 
odically travel past our line of sight as the star rotates. The 
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pulsation period is equated with the rotation period of the 
neutron star’s crust and this value changes with time as ac- 
creted matter transfers angular momentum to or from the 
star. 

HMXB are composed from massive, and therefore rel- 
atively short lived, stars and so their presence indicates 
that a new population of stars formed relatively recently 
(Grimm, Gilfanov and Sunyaev 2003). This may be why 
more HMXB have been discovered in the Milky Way and 
the Small Magellanic Cloud (SMC) than the Large Magel- 
lanic Cloud (LMC). 

The HMXB in the Milky Way are mostly found 
within the spiral arms (Coleiro and Chaty 2011) and the 
SMC underwent a burst of star formation about 200 mil- 
lion years ago when it was about three times closer to 
the LMC than it is today (Gardiner and Noguchi 1996; 
McBride et al 2010). The SMC has an HMXB popula- 
tion which is comparable in number to the Galaxy de- 
spite being almost two hundred times less massive. Whilst 
about half of the HMXB in the Galaxy are BeXB, all 
but one HMXB in the SMC is a BeXB system. The 
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Observation ID 

Start Date 

Start Time 

Exposure (ks) 

00038029001 

2008-10-23 

03:56:38 

6.56 

00038029002 

2008-11-11 

14:35:46 

6.40 

00038029003 

2008-11-14 

00:18:15 

2.59 

00038029004 

2009-01-10 

13:23:59 

4.32 

00038029005 

2011-09-30 

16:50:35 

4.47 

00038029006 

2011-10-03 

18:28:46 

5.90 


Table 1. Summary of XRT datasets. 


LMC, on the other hand, contains relatively few HMXB, 
most of which are BeXBs (Negueruela and Coe 2002; 
Liu, van Paradijs and van den Hauvel 2005; Sturm et al 
2012 ). 

Swift J045106. 8-694803 is part of a newly discovered 
BeXB system located in the LMC. It was first discov- 
ered by Beardmore et al (2009) at an RA, Dec (J2000) 
of 04:51:06.8 and -69:48:03.2 respectively, with an uncer- 
tainty of 3. 5". Beardmore et al measured an X-ray flux of 
(1.68±0.11)xl0‘ 11 erg cm' 1 2 s 4 over 0.3-10 keV, fit with a 
power law of photon index 0.96±o;ot with a column density 
of (1.9±0.3)xl0 21 cm' 2 . 

The 14-195 keV flux was found to be (2.8±0.3) x 10' 11 
erg cm' 2 s _1 fit with a power law of photon index 2.5±0.4. 
Beardmore et al also report a pulsation period of 187 sec- 
onds. The optical companion to this X-ray source is the V = 
14.70 blue star [M2002] 9775 (Massey 2002) with an orbital 
period of 21.64T0.02 days. 

Grebenev et al have recently calculated a photon index 
of 0.5T0.5 using data from INTEGRAL. They have also cre- 
ated energy spectra over 3-200 keV showing that the energy 
of the high energy cut-off in the spectrum is at 16.0T5.0 keV 
(Grebenev, Lutovinov and Tsygankov 2012). 

In this paper we look at Swift J045106. 8-694803 in more 
detail, confirming its luminosity, pulse period, orbital pe- 
riod and spin-up rate which can be used to determine the 
strength of its magnetic field. 

2 OBSERVATIONS 

2.1 X-ray Observations 

2.1.1 Swift/ XRT 

Swift’s X-ray telescope (XRT) is a CCD imaging spectrom- 
eter operating over 0.2-10 keV in photon counting mode. 
Archival data were taken from NASA’s High Energy As- 
trophysics Science Archive Research Center (HEASARC) 1 
covering six time periods from 23rd October 2008 to 3rd 
October 2011, as summarised in Table 1. 

The images were extracted using the ftool 2 xselect. 
Source and background spectra were then extracted from 
regions of 34 ,, radii. The spectra were binned to have 50 
counts per bin. The ancillary response files (ARF) were 
calculated with xrtmkarf and a redistribution matrix file 
(RMF) was taken from HEASARC’s calibration database 
(CALDB). The position of the source was confirmed using 
ftool xrtcentroid. 

The total count rate and error of each dataset, as 

1 http://heasarc.gsfc.nasa.gov/ 

2 http: / /heasarc. nasa.gov/ftools/ 


Observation ID 

Start Date 

Start Time 

Exposure (ks) 

96441-01-01-00 

96441-01-01-01 

2011-10-28 

2011-10-28 

13:54:55.6 

10:51:38.6 

9.005 

0.849 


Table 2. Summary of RXTE datasets. 

well as the intrinsic H I column density (Nh) and pho- 
ton index, were calculated using ftool xspec. The spec- 
tra were described by an absorbed power law with a 
fixed Galactic foreground column density of 8.4 xlO 20 cm" 2 
(Dickey and Lockman 1990) and abundances set in accor- 
dance with Wilms, Allen and McCray (2000). Intrinsic ab- 
sorption and the abundances of elements heavier than he- 
lium were set to 0.4 (Borkowski, Hendrick and Reynolds 
2006). X-ray spectra were then compiled in xspec over 0.2-10 
keV. 

The 0.2-10 keV flux of each dataset was determined 
using xspec. The luminosity was then calculated using a dis- 
tance of 50.6±1.6 kpc (Bonanos et al 2011). 

The light curve of each dataset was extracted in xselect 
and a Lomb-Scargle normalised periodogram was produced 
for each using time-series analysis package Period. 3 . 

2.1.2 Swift /BAT 

The Swift Burst Alert Telescope (BAT) data were taken 
with the Hard X-ray Transient Monitor from 16th December 
2004 to the 31st May 2010 over 14-195 keV. A 58 month 
light curve was downloaded from NASA’s Swift/BAT 58- 
Month Hard X-ray Survey 4 . This contained an average of 
~15 observations a day split into 8 energy bands. 

2.1.3 RXTE 

Archival data were taken from HEASARC. These were 
recorded in two datasets on the 28th October 2011 using 
the Rossi X-ray Timing Explorer’s (RXTE) Proportional 
Counter Array, over ~3-10 keV, as summarised in Table 2. 

Cleaned light curves were produced for each dataset 
and combined. A Lomb-Scargle normalised periodogram was 
then created using a frequency interval of lxlO" 5 Hz. 

2.2 Optical observations 

2.2.1 Optical photometry 

Archival data were taken from the MAssive Compact Halo 
Objects (MACHO) project using the 1.27 meter telescope at 
the Mount Stromlo Observatory in Australia. This covered 
the period from 1st November 1992 to the 29th December 
1999 and contains instrumental magnitudes using red (R) 
and blue (B) filters. 

The data were filtered to remove results flagged as er- 
roneous. The four points remaining in the R. band dataset 
which were over 2.4 standard deviations from the mean were 
also removed. This left 96 data points in the red and 207 in 
the blue. A Lomb-Scargle normalised periodogram was then 
created in Period for both the R and B band datasets using 

3 http:/ /www. starlink.rl.ac.uk/docs/sunl67.htx/sunl67.html 

4 http:/ /swift. gsfc.nasa.gov/docs/swift/results/bs58mon/ 
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a frequency interval of lxlO' 5 Hz and detrending using a 
polynomial of order 3. 


2.2.2 Optical spectroscopy 

Optical spectra of Swift J045106. 8-694803 have been taken 
on three separate occasions. Red end spectra were taken 
with the 1.9m Radcliffe telescope at the South African As- 
tronomical Observatory (SAAO) on the 12th December 2009 
and 26th September 2011. The data were obtained using the 
unit spectrograph combined with a 1200 lines mm” 1 grating 
and the SITe detector at the Cassegrain focus. The resulting 
spectra have a spectral resolution of ~ 3 A. Comparison cop- 
per neon spectra were taken immediately before and after 
the observation and were used for wavelength calibration. 

The source was also observed on the 8th and 10th De- 
cember 2011 with the New Technology Telescope (NTT), 
La Silla, Chile. Grisms 14 and 20 on the ESO Faint Ob- 
ject Spectrograph (EFOSC2) were used for blue and red 
end spectroscopy respectively with a slit width of 1.5”. 
Grism 14 has a grating of 600 lines mm -1 that yielded 
1 A pixel -1 dispersion over a wavelength range of AA3095- 
5085 A. Grism 20 is one of the two new Volume-Phase Holo- 
graphic grisms recently added to EFOSC2. It has a smaller 
wavelength range, from 6047-7147 A, but a superior disper- 
sion of 0.55 A pixel -1 and 1070 lines pixel -1 . Filter OG530 
was used to block second order effects. The resulting spectra 
were recorded on a Loral/Lesser, thinned, AR-coated, UV 
flooded, MPP CCD with 2048x2048 pixels, at a spectral 
resolution of ~ 10 A and ~ 6 A respectively. Wavelength 
calibration was achieved using comparison spectra of He- 
lium and Argon lamps taken through out the observing run 
with the same instrument configuration. 

The data were reduced using the standard packages 
available in the Image Reduction and Analysis Facility 
(IRAF). The resulting spectra were normalized to remove 
the continuum and and a redshift correction applied corre- 
sponding to the recession velocity of the LMC, taken as -280 
km s -1 (Paturel et al 2002). 



Figure 1. Uncertainty circles of the X-ray source overlaid on an 
optical background (from ESO DSS II Blue). The red circle shows 
the position found using the Swift/XRT dataset 00038029001, the 
blue circle using dataset 00038029002 and the green using the 
XMM-Newton dataset from Bartlett et al (2012, in prep). 



Figure 2. The upper panel shows the X-ray spectrum from 
Swift/XRT dataset 00038029001 over 0.2-10 keV, here the pho- 
ton index is 0.7±0.1 with a Njj of 3±2 counts/s, the lower panel 
shows residuals. 


3 RESULTS 

3.1 X-ray Results 

3.1.1 Swift/XRT 

Figure 1 shows the positions calculated for datasets 
00038029001 (red) and 2 (blue). The first dataset has an RA, 
Dec (J2000) of 04:51:06.4 and -69:48:02.5 with a two sigma 
error radius of 3.6” . The second has an RA, Dec (J2000) of 
04:51:07.0 and -69:48:03.1 with a two sigma error radius of 
3.6”. These are consistent with the positions calculated by 
Beardmore et al (2009). 

The results of data extracted using xspec are sum- 
marised in Table 3. Figure 2 shows the X-ray spectra of 
dataset 00038029001 over 0.2-10 keV. This shows that the 
counts can be divided equally into four energy ranges of 0.5- 
1.5, 1.5-3, 3-4.5 and 4.5-8 keV. The X-ray luminosities are 
plotted in Figure 3, which shows that the X-ray luminosi- 
ties measured between October 2008 and January 2009 were 
about four times higher than those measured in October and 
September 2011. 


Figure 4 shows the Lomb-Scargle normalised peri- 
odogram for each of the six datasets. A period is only ev- 
ident in datasets 00038029001 and 2, these are 187.0T0.3 
and 186.8±0.2 seconds respectively, and show no harmonics. 
This indicates that light curves folded at these values should 
be sinusoidal. Folded light curves for datasets 00038029001-2 
are shown in Figure 5. 

Light curves for dataset 00038029001 in each of the four 
energy ranges determined from Figure 2 are shown in Figure 

6. A plot of pulse- fraction against energy is shown in Figure 

7. The pulse-fraction appears to increase with increasing en- 
ergy between 0.5 and 4.5 keV and then decrease, however 
the large error bars make this inconclusive. 

The side peaks in the periodograms are due to gaps in 
the observations. This was confirmed firstly, by fitting a sine 
wave with the same period to the data, which gave rise to 
identical peaks. Secondly, by splitting dataset 00038029001 
into 6 shorter datasets composed of continuous observations 
and adding the individual periodograms, as can be seen in 
Figure 8. 
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Observation ID 

Count Rate 
(counts / s 
over 0.2-10 keV) 

N h 

(10 21 cm" 2 ) 

Photon 

Index 

Flux 

(10 -11 ergs/ 
s/cm 2 ) 

Luminosity 
(10 36 ergs/s) 

Reduced 

Chi-Squared 

Degrees of 
freedom 

a 

00038029001 

0.1984=0.006 

3±2 

0.74=0.1 

2.34=0.1 

7.04=0.3 

1.4 

21 

b 

00038029002 

0.1754=0.005 

1±2 

0.74=0.1 

2.14=0.1 

6.34=0.4 

1.2 

17 

c 

00038029003 

0.1824=0.008 

3±7 

0.64=0.2 

2.34=0.2 

7.04=0.6 

0.6 

5 

d 

00038029004 

0.1814=0.006 

1±3 

0.64=0.1 

2.14=0.1 

6.54=0.4 

1.2 

11 

e 

00038029005 

0.0734=0.004 

15±14 

1.14=0.4 

0.74=0.2 

2.24=0.5 

0.8 

2 

f 

00038029006 

0.058±0.003 

14=16 

0.64=0.4 

0.84=0.2 

2.34=0.5 

2.0 

2 


Table 3. Summary of information extracted from the spectra of Swift/XRT datasets 00038029001-6 over 0.2-10 keV using a distance of 
50.6d=1.6 kpc (Bonanos et al 2011). 
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Figure 3. Luminosity plotted against time for the Swift/XRT 
and XMM-Newton datasets over 0.2-10 keV. XMM-Newton 
dataset from Bartlett et al (2012, in prep). 


Figure 5. Light curves from Swift/XRT datasets 00038029001 
(top panel) and 00038029002 (bottom panel) over 0.2-10 keV, 
folded at 187.0 and 186.8 seconds respectively. The phase shift is 
arbitrary. 



Peril :l ! sees ) P; riocl (s.ecs) 


Figure 4. Lomb-Scargle normalised periodograms for Swift/XRT 
datasets 00038029001-6 (labeled a-f) over 0.2-10 keV. The dotted 
lines indicate the 187 second period and the two side peaks of 181 
and 193 seconds mentioned by Beardmore et al (2009). 

Datasets 00038029001 and 2 provided the best signal 
and so their light curves were combined to determine the 
best estimation for the pulse period in late 2008. A Lomb- 
Scargle normalised periodogram yielded a maximum peak 
at 186.85±0.04 seconds as can be seen in Figure 9. Monte 
Carlo simulations give this period a 99.9% confidence level. 




as 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0 

PI ase Si ise 

Figure 6. Light curves from Swift/XRT dataset 00038029001, 
folded at 187.0 seconds, over 0.5-1. 5, 1.5-3, 3-4.5 and 4.5-8 keV 
(labelled a-d). 


Bootstrapping was conducted in order to confirm that this 
is the correct peak, and not a product of the window func- 
tion, and after 5000 iterations between 166 and 200 seconds, 
a period of 187.074=0.04 seconds was found 63% of the time 
and 181.41±0.04 seconds 37% of the time. We therefore con- 
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Figure 7. Pulse-fraction plotted against energy for Swift/XRT 
dataset 00038029001 folded at 187.0 seconds, over 0.5-1. 5, 1.5-3, 
3-4.5 and 4.5-8 keV. 
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Figure 8. Lomb-Scargle normalised periodogram created by 
splitting dataset 00038029001 into 6 datasets composed of contin- 
uous observations and adding the periodograms. The dotted line 
indicates the 187.0 second period found for dataset 00038029001. 

eluded that 187.07±0.04 seconds was the most probable pe- 
riod of Swift J045106. 8-694803 in October/November 2008. 

3.1.2 Swift/B AT and INTEGRAL/IBIS 

The top panel of Figure 10 shows the long-term light curve of 
Swift J045106. 8-694803 with total counts over 14-195 keV, 
binned at 28 days. The count rate appears to increase over 
time, peaking in about July 2007 and continuing to remain 
above zero for the next three years. 

This was confirmed by the light curve compiled from the 
INTEGRAL IBIS data which is shown in the second panel 
of Figure 10. The bottom panel shows the light curve of the 
optical counterpart taken from OGLE III and IV in the I- 
band. The OGLE data shows that the source was brightest 
in the I-band whilst it was barely detectable in the X-ray. 
There is a slight increase in brightness after MJD 54000 but 
this is not significant. 

A Lomb-Scargle normalised periodogram was created 


CD 

i 

o 

CL 


Figure 9. Lomb-Scargle normalised periodogram from 
Swift/XRT datasets 00038029001 and 2 over 0.2-10 keV 
showing a maximum peak at 186.85i0.04 seconds. The side 
peaks are due to gaps in the observations. 


for the BAT dataset using Period with a frequency interval 
of 1 x 10' 5 Hz but no periods were detected. 


3.1.3 RXTE 

The Lomb-Scargle normalised periodogram is shown in plots 
(a) and (b) of Figure 11. As with the Swift/XRT data, the 
side peaks are due to gaps in the observations. This was 
confirmed by fitting a sine wave with the same period to the 
data, as can be seen in plot (c) of Figure 11. The harmonics 
at exactly 1/3 and 1/4 of the pulse period, seen in plot 
(a), indicate that the light curve should be non-sinusoidal 
and have multiple peaks when folded at the pulse period, as 
shown in Figure 12. 

The shape of the light curve contains information on 
the emission geometry from the regions close to the neutron 
stars magnetic poles. At high luminosities photons can es- 
cape from the sides of the accretion column, giving rise to 
a fan-beam pattern and a more complex profile than that 
found at lower luminosities where X-rays are generally emit- 
ted in a pencil-beam. 

The RXTE data show a pulsation period of 169.8±0.3 
seconds, 17.3±0.3 seconds less than the period evident in 
the Swift/XRT datasets taken in October and November 
2008. This gives a spin-up rate ( P ) of -5.8±0.1 s/yr, or 
(-I.84i0.03) x 10' 7 s/s, and a spin-up time scale (T s ) of - 
30.8A0.5 yr using 



where P is the average pulse period in seconds. Since this is 
an average over three years, it is possible that the spin-up 
rate has been much higher at points during this period. 
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Figure 10. Long-term light curve from Swift/BAT (14-195 keV), INTEGRAL IBIS (15 keV - 10 MeV) and OGLE III and IV. 



Pedod (secs) F'eried (secs) 


Figure 11. Panel (a) shows the Lomb-Scargle normalised peri- 
odogram from the combined RXTE datasets over 3-10 keV. Panel 
(b) shows a close up view of the region around the main peak com- 
pared to panel (c) which shows the result of a simulated dataset 
produced by a pure sine wave of period 169.78 seconds. 


3.1.4 XMM-Newton 

Swift J045106. 8-694803 was observed by the X-ray Multi- 
Mirror Mission - Newton (XMM-Newton) on 17th July 2012 
(Bartlett et al 2012, in prep). This confirmed the position of 
Swift J045106. 8-694803 - as is shown in Figure 1 - and gave 



Phase 


Figure 12. Light curve folded at 169.78 seconds from RXTE 
datasets 96441-01-01-00 and 01 over 3-10 keV. 


an up-to-date period and luminosity. The position was found 
to be at an RA, Dec (J2000) of 04:51:06.7 and -69:48:04.2 
respectively, with a one sigma uncertainty of 1". The lumi- 
nosity was found to be (9.8±0.9) x 10 34 ergs/s - as shown in 
Figure 3 - and the period 168.5±0.2 seconds. 

This is 1.3±0.4 seconds less than the period calculated 
from RXTE in 2011, giving a P of -1.8±0.5 s/yr. It is 
also 18.6±0.2 seconds less than the period calculated from 
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Figure 13. Plot showing the three periods measured by 
Swift/XRT in 2008, RXTE in 2011 and XMM-Newton in 2012 
(Bartlett et al , 2012, in prep). The straight line is a line of best 
fit showing that the period is continuing to decrease, despite the 
fact that the rate of change is also decreasing. 


Figure 15. Lomb-Scargle normalised periodogram for the B 
band MACHO dataset showing a possible orbital period of 
21.631±0.005 days. 
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Figure 14. Light curves for the R and B band MACHO datasets. 
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Figure 16. Light curve from the B band MACHO dataset folded 
at 21.631 days. 


Swift/XRT in 2008, giving an average P of -5.01±0.06 s/yr 
and a T s of -34.9±0.4 yr. Figure 13 shows the three periods 
measured by Swift/XRT in 2008, RXTE in 2011 and XMM- 
Newton in 2012. The luminosity and rate of change of spin 
period may be decreasing but Swift J045106. 8-694803 is still 
continuing to spin-up at a high rate. 

3.2 Optical results 

3.2.1 Optical photometry 

Figure 14 shows that the flux of the optical companion to 
Swift J045106. 8-694803 appears to have remained fairly con- 
sistent in the B and R bands over 7 years. The B band 
dataset show a possible orbital period of 21.631±0.005 days 
as seen in Figure 15, although any underlying non-radial pul- 
sations from the Be star may affect the results (Bird et al 
2012). Non-radial pulsations occur when some parts of the 
stellar surface move inwards while others move outwards at 
the same time. It has been suggested that they could help 


form the Be stars circumstellar disk. If the rotational ve- 
locity of the Be star is close to the critical velocity, then 
pulsations lead to matter being ejected and spun up to form 
a Keplerian disk (Reig 2011). 

This period is not evident in the R band, this is most 
likely due to the lack of data points rather than the conse- 
quence of a real effect. This idea was confirmed by randomly 
removing half of the B band data points and creating a new 
Lomb-Scargle normalised periodogram, which also failed to 
show any evidence of an orbital period. A light curve folded 
at 21.631 days was then produced for the B band dataset as 
shown in Figure 16. Colour ratio and colour magnitude dia- 
grams were also produced. These confirmed that the magni- 
tude in each colour band has remained fairly consistent over 
the observation period. 

Swift J045106. 8-694803 was also observed as part 
of the OGLE III (LMC136.6.14874) and OGLE IV 
(LMC531.05.4251) programmes. These I-band data cover, 
in total, a duration of more than a decade, and are shown 
in the bottom panel of Figure 10. From this figure it is clear 
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that the source undergoes a significant long-term modula- 
tion on periods in excess of 400 days. Formally, a Lomb- 
Scargle power spectrum gives the peak to be around 440 
days, but a visual inspection of the light curve shows that 
there are other time-scale changes occurring. It is therefore 
unlikely that this long period is related directly to the binary 
period, but rather may either indicate a general time-scale 
for fluctuations in the stellar wind, or precessional motion 
of the circumstellar disk. All these variations make it very 
difficult to search for confirmation of the 21.631 day pe- 
riod seen in the MACHO data. However, if only the better- 
sampled OGLE IV data are merged with the MACHO data 
(normalised to the approximate starting magnitude of the 
OGLE III data), then the strength of the 21.631 day peak 
in the Lomb-Scargle power spectrum increases slightly. But 
without the prior knowledge of this period from the MA- 
CHO data such a period would not have been found in the 
OGLE data. 

3.2.2 Optical spectroscopy 

OB stars in the Milky Way are classified using certain metal 
and helium line ratios (Walborn and Fitzpatrick 1990) 
based on the Morgan-Keenan system (MK; Morgan et al 
1943). However, this is unsuitable in lower metallicity en- 
vironments as the metal lines are either much weaker or not 
present. As such, the optical spectrum of IGR J05414-6858 
was classified using the method developed by Lennon (1997) 
for B-type stars in the SMC and implemented for the SMC, 
LMC and Galaxy by Evans et al (2004, 2006). 

Figure 17 shows the unsmoothed optical spectrum of 
Swift J045106. 8-694803. The spectrum is dominated by the 
hydrogen Balmer series and neutral helium lines. There does 
appear to be evidence for the Hell A4200A line, but it is dif- 
ficult to distinguish above the noise level along with the Hell 
AA4541, 4686A lines. The Hel A4143A line is clearly stronger 
than the Hell A4200A line constraining the optical counter- 
part of Swift J045106. 8-694803 to be later than type 09. 
There also appears to be evidence for the Siiv A4088, A4116A 
lines necessary for a B1 classification. This is supported by 
the relative strengths of the Sim A4553A and Mgn A448lA 
suggesting a classification of B2 or earlier. 

The luminosity class of the companion star was de- 
termined by the ratios of SlV A4088/Hel AA4026 — 412lA, 
Siv A4116/Hei A4121 and Hen A4686/Hei A4713. The first 
two ratios strengthen with decreasing luminosity class (i.e. 
with increasing luminosity) whereas the latter ratio de- 
creases with increasing luminosity. The relative strengths 
of these lines are contradictory: The Hell A4686/Hel A4713 
and SlV A4116/Hel A4121 ratios suggest a luminosity class 
III, although the proximity of the rotationally broadened 
H<5 A4102A to the Siiv lines makes this more complex. The 
Siv A4088/Hei A4026 ratio is more consistent with a star 
of luminosity class V. The V band magnitude of this star 
is reported by several sources as between 14.6 and 14.7 
(e.g. Massey 2002; Zaritsky et al 2004). An My of -4.3±0.1 
was calculated using a distance modulus of 18.52±0.07 
(Bonanos et al 2011) along with an my of 14.65±0.05 and 
an extinction Av=0.4±0.1 (calculated using the Galactic 
column density towards the source, 8.4 ± 1.0 x 10 20 cm ” 2 
and the results of Giiver and Ozel 2009). This is consistent 
with B0.5III (Wegner 2006). Less information is available for 


the absolute magnitude of emission line stars, which will be 
dependent on the inclination and size of the disc. As such we 
classify the optical counterpart of Swift J045106. 8-694803 as 
a B0-1 III-V star. 

Figure 18 shows the three red end spectra of Swift 
J045106. 8-694803 taken two years apart. The ESO (top) 
spectrum is offset by 6 flux units. The Ha equivalent width, 
considered an indicator for circumstellar disk size, is remark- 
ably similar in all three spectra increasing from -29±2 A 
for the SAAO spectrum taken in 2009 to -33±1 A and - 
34.5±0.6 A for the SAAO and ESO spectra taken in 2011. 
It is not uncommon to see large variations in the Ha equiv- 
alent width of these systems on timescales of months. This 
apparent consistency in the equivalent width (and hence the 
disk size) is almost certainly linked to the exceptionally per- 
sistent X-ray activity of the source. 


4 DISCUSSION AND CONCLUSIONS 

We have confirmed the pulsation and orbital periods given 
by Beardmore et al (2009) . Recent measurements show that 
Swift J045106. 8-694803 has spun-up at a rate of ~5 seconds 
a year for the last four years and although this rate is de- 
creasing with luminosity, it is still continuing to spin-up at 
a high rate. 

Swift J045106. 8-694803 may have experienced one of 
the highest spin-up rates of any known accreting pulsar but 
this can be accounted for using Ghosh and Lamb’s accretion 
theory (1979) assuming that the stellar magnetic field is 
dipolar and that it is accreting from a Keplerian disk. 

Accretion theory shows that the magnetic field (B) of 
a disc-fed accreting pulsar is related to the rate of change 
of its spin period (P), its average spin period (P) and its 
average luminosity (L). 

In general, the higher the mass accretion rate ( M ), the 
higher the luminosity and the more angular momentum is 
accreted per second. Where magnetic stresses dominate mat- 
ter flow at surface SI - as seen in Figure 19 from Ghosh and 
Lamb (1979) - this almost always causes spin-up. 

At surface S2, where viscous stresses dominate, this can 
add a spin-up or spin-down torque depending on the stars 
fastness parameter (<u s ). This is the ratio between the angu- 
lar velocity of the neutron star and the Keplerian angular 
velocity of the disc. 

If the pulse period is long, plasma is able to penetrate 
into the magnetosphere and reach the neutron star surface, 
transferring angular momentum and causing a strong spin- 
up torque. If the pulse period is short, the plasma is unable 
to penetrate further and the neutron star spins down. 

The magnetic moment (p) can be found if P, P and L 
are known by assuming a given mass and radius and seeing 
which value of fi predicts the observed P using 

- P = 5.0 x 10~ 5 l 4 / 0 7 n(u s )S 1 (M)(PL 3 3 / 7 7 ) 2 ( 2 ) 

where p, 3 o is the magnetic moment in units of 10 30 Gauss 
cm ' 3 and P 37 is the luminosity of the accreting star in units 
of 10 37 ergs/s. For an 0< w s <0.9 

n(u] s ) = 1.4(1 — 2.86tu s )(l — w ,,)” 1 (3) 
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Figure 17. Spectrum of Swift J045106. 8-694803 in the wavelength range AA3900-5000A with the NTT on 2011-12-08. The spectrum has 
been normalized to remove the continuum and redshift corrected by -280 km s — 1 . Atomic transitions relevant to spectral classification 
have been marked. 



Wavelength (A) 

Figure 18. ESO (top panel) and SAAO (middle and bottom 
panels) spectra of Swift J045106. 8-694803 in the wavelength range 
AA6400-6700A with the NTT on 2011-12-10. Spectra have been 
normalized to remove the continuum and shifted by -280 km s _1 


within 5% accuracy and 

UJ S = 1.35^% 7 S 2 (M)(PLl / 7 7 )- 1 = 1.19 

( 4 ) 


where M17 is the mass accretion rate in units of 10 1 ' grams/s 
and is assumed throughout to be 1.4. S'i(M) and S 2 (M) 
are structure functions that depends on the mass equation 
of state and the dynamical response of the neutron star. 

Si(M) = Re 7 (^)~ 3/7 I^ (5) 

where Re is the radius of the accreting star in units of 10 6 cm, 
assumed throughout to be 1, and 745 is the moment of inertia 
in units of 10 45 grams cm 2 . 

S 2 (M) = T?- 3/7 (^-)- 2/7 (6) 

The value of the magnetic field can then be calculated using 


Since n(w a ) depends only on w a and w s is inversely pro- 
portional to PL37 3 / 7 , the P for a star of a given mass and 
magnetic moment depends only on (PL37 3 / 7 )" 1 . Figure 20 
shows that Swift J045106. 8-694803 is located in the same 
region as GX 1+4, SAX J2103. 5+4545, 4U 2206+54 and 
SXP 1062 on a plot of logio(.P) against logio(PL37 3 ^ 7 ). 

GX 1+4 is part of a LMXB, it is a slow rotator 
with a period of about 121 seconds and has been ob- 
served to have spun-down by about 2.6 seconds a year 
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Figure 19. Side view of accretion flow from Ghosh and Lamb (1979). 


(Chakrabarty et al 1997). SAX J2103. 5+4545 is part of a 
HMXB system, it is also a slow rotator with a pulse period of 
358.166+0.0005 seconds and a luminosity of (2.0+0.5) x 10 36 
ergs/s. It has also been known to spin- up at a rate of 
4.69+0.09 seconds per year (Sidoli et al 2005). This means 
that it has a magnetic field of (1.47±o;og) x 10 14 Gauss. 
4U 2206+54 (Reig, Torrejon and Blay 2012) and SXP 1062 
(Turolla and Popov 2012) are also part of HMXB systems, 
both are known to be spinning down and both were recently 
identified as having extremely high magnetic fields. 

Highly magnetised neutron stars may be formed if a 
very magnetic star conserves its magnetic field during the 
neutron stars formation or, if the neutron star is rotating 
fast enough to start a dynamo effect during the first ~20 
seconds of its life. This converts heat and rotational en- 
ergy into magnetic energy and increases the magnetic field 
(Ferrario and Wickramasinghe 2006). 

The first evidence for highly magnetised accreting neu- 
tron stars came from Pizzolato et al (2008), who showed 
that the X-ray source 1E161348-5055 may contain a neu- 
tron star with a surface magnetic field of ~10 15 G which is 
part of a LMXB. This was shortly followed by Bozzo et al 
(2008), who argued that highly magnetised neutron stars 
may also exist in HMXB. 

Our results show that the P measured between 2008 
and 2011 implies a magnetic field of (1.3+0. 1) x 10 14 Gauss, 
whilst that measured between 2011 and 2012 implies a mag- 
netic field of (6+1) xlO 13 Gauss. The value calculated for 
the whole duration of observations is (1.2+0. 1) x 10 14 Gauss. 
Taking into account the highest and lowest possible values of 
B measured between 2008 and 2011 and 2011 and 2012, this 
value becomes (1.2+0;!) x 10 14 Gauss. Parameters calculated 
using equations 1-7 are shown in Table 4. 

All values are over the quantum critical value of 
4.4xl0 13 Gauss, which means that the physics must be 
described with quantum field theory rather than classical 
physics. Although Swift J045106. 8-694803 is powered by ac- 
cretion, rather than due to the decay and instabilities of its 
magnetic field, and hence cannot be referred to as a mag- 
netar yet, it has the potential to become one. The Galactic 


Orbital period 
Pulse period in October 
/November 2008 
Pulse period in October 2011 
Pulse period in July 2012 
Average Luminosity 

21.631+0.005 days 
187.07+0.04 s 

169.8+0.3 s 
168.5+0.2 s 
(3.4±0.3)xl0 36 ergs/s 

p 

-5.01+0.06 s/yr 
(-1.59+0.02) xlO' 7 s/s 

T s 

-34.9+0.4 yr 

logic (P) 

0.70+0.01 

logic (PL 37 3 / 7 ) 

2.04+0.04 

UJs 

0-7+^ 

n(u) s ) 

- 3 ±i 

h 

(1.2±q !)x 10 32 Gauss cm' 3 

B 

(1.2±° !)xl0 14 Gauss 


Table 4. A summary of observed data and results derived using 
equations 1-7. for Swift J045106. 8-694803. 


HMXB system containing LSI 61303 recently showed mag- 
netar like behaviour when it underwent two bursts similar 
to those of soft gamma repeaters (SGRs) (Torres et al 2012; 
Papitto, Torres and Rea 2012). 

In conclusion, the spin-up rate of Swift J045106.8- 
694803 can be explained using Ghosh and Lamb’s (1979) 
accretion theory provided that it has a magnetic field of 
(1.2 ±o;!)x 10 14 Gauss. 
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Figure 20. Logarithmic plot of P against PLr with contours calculated from Dachs, Hummel and Hanuschik (1992). Left panel shows 
the theoretical curves for three values of stellar magnetic moment assuming a neutron star mass of 1.3 Mg. The dashed line shows 
the theoretical curve for //30 = 0.48. The right panel shows the theoretical curves for three values of neutron star mass assuming a 
magnetic moment of //.3Q = 0.48. The dashed line is the theoretical curve for M=1.3 Mg. Data for Swift J045106. 8-694803 and seven 
other pulsars are included in the plot. References are as follows; Her X-l (Li, Dai and Wang 1995), SMC X-l (Kahabka and Li 1999), X 
Per (Delgado-Marti et al 2001), GX 1+4 (Ferrigno et al 2007 and Chakrabarty and Roche 1997), SAXJ2103. 5+4545 (Sidoli et al 2005), 
4U 2206+54 (Reig, Torrejon and Blay 2012) and SXP 1062 (Henault-Brunet et al 2012 and Haberl et al 2012). 
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